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Summary. We report here on angular broadening measurements of 20 
OH/IR stars near the Galactic center. This class of sources is known to 
have bright, intrinsically compact (< 20 mas) maser components within 
their circumstellar shells. We have used VLB A antennas and the VLA 
to perform a MKII spectral line VLBI experiment. The rapid drop in 
correlated flux with increasing baseline, especially for sources closest to 
the Galactic center, is attributed to interstellar scattering. Angular di- 
ameters were measured for 13 of our sources. Lower limits were obtained 
for the remaining 7. With our data, together with additional data taken 
from the literature, we have been able to determine the distribution of 
interstellar scattering towards the Galactic center. We find a region of 
pronounced scattering ( 6 <*, = 525 mas at 1.612 GHz) nearly centered 
(angularly) on SgrA*. We consider two interpretations for the enhanced 
scattering. One hypothesis is that the scattering is due to a clump of 
enhanced turbulence, such as those that lie along lines of sight to other 
known objects, that has no physical relationship to the Galactic Cen- 
ter. The other model considers the location of the enhanced scattering 
to arise in the galactic center itself. The physical implications of the 
models yield information on the nature of interstellar scattering. Future 
tests to discriminate between these models will be the measurement of 
angular sizes of extragalactic sources within 15' of SgrA* and angular 
sizes of additional OH/IR stars toward the Galactic center region. 
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1. Introduction 

It has long been suspected that the observed size of the Galactic center source SgrA* 
is set by interstellar scattering (Davies et al., 1976, Lo et al., 1981, Lo et al., 1985, 
Backer, 1988). At 1.6 GHz, the frequency we use throughout this paper, the scattered 
broadened diameter of SgrA* is about 520 mas. Such a high value of interstellar 
scattering is unusual, as there are only a few lines of sight with comparable degrees of 
plasma turbulence (Cordes, Ananthakrishnan, and Dennison 1984; Moran et al., 1990, 
Fey et al., 1991). Is this scattering anomalous, local to SgrA* and a consequence of the 
activity occurring there (Backer, 1988)? Is it due to clumps of strong turbulence at 
some random position along the line of sight (Cordes, Weisberg, and Boriakolf, 1985)? 
Or can it be explained as part of the general increase in scattering that is observed 
towards the inner Galaxy (Fey et al., 1991, Rao & Slee, 1988)? The answer to this 
question has a bearing on an unsolved problem in the study of plasma turbulence in 
the interstellar medium: the identity of the agent(s) responsible for the turbulence. 

Recently, van Langevelde & Diamond (1991) showed that the OH masers in the 
circumstellar shells of four OH/IR stars, all less than 16 arcmin from SgrA*, were 
broadened by an amount similar to SgrA*. Thus they confirmed that the size of 
SgrA* at centimeter wavelengths is set by interstellar scattering. An equally impor- 
tant consequence of this work was the realization that OH/IR stars are a powerful 
tool to study the turbulent electron density component in the very inner regions of 
our Galaxy. The number of observable OH/IR stars in this region is at least an order 
of magnitude larger than that of extragalactic sources, and they are all at roughly 
similar distances. Moreover, as we will show in section 2, there exist bright, intrinsi- 
cally compact (< 20 mas) maser spots in all OH/IR stars and so the usual problem 
of intrinsic source structure for extragalactic sources does not enter into the inter- 
pretation of the visibility data. Pulsars, which in general are excellent point-source 
probes of interstellar scattering, cannot be used to obtain information on interstellar 
scattering close to the Galactic center. Although pulsars are probably present within 
a few kpc of the Galactic center (e.g. Clifton & Lyne 1986), it is unlikely that they 
will be found within 100 pc, because interstellar scattering will quench the pulsations. 

In this paper we expand on the earlier work of van Langevelde ic Diamond (1991) 
with the goal of determining the distribution of the scattering material around SgrA* 
in order to place constraints on the properties and the location of this medium. The 
outline of this paper is as follows. We start by describing our observational technique 
as it applies to the angular broadening of OH/IR stars specifically. This includes 
setting up the necessary formulas to describe scattering as well as the relation between 
interstellar scattering and the emission measure, which we will use to discuss the 
possibility of a scattering medium located at the Galactic center. In section 3 we 
discuss our source selection criteria, the observational parameters, and the reduction 
of the visibility data to scatter diameters. Finally, in section 4, we use these diameters, 
in addition to those from the literature, to examine possible models for the source of 
scattering in the direction of the Galactic center. 
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2. Technique 

2.1. ANGULAR BROADENING BY INTERSTELLAR SCATTERING 

Angular broadening has its origin in fluctuations of the electron density in the in- 
terstellar medium. Because the refraction index of the interstellar medium for radio 
waves is given by n « 1 - (r e n e A 2 /2jr), fluctuations in n a cause disturbances in a 
wave front travelling through the interstellar medium. This degrades the degree of 
coherence between the signals received at two widely separated points, by an amount 
that depends on the level of turbulence present in the medium. Electron density fluc- 
tuations on scales of 10® - 10® m are readily measurable with an interferometer by 
determining either the loss of coherence as a function of antenna separation (in the 
visibility plane) or, equivalently, the angular broadening (in the image plane). 

The same fluctuations give rise to a variety of other observational phenomena (see 
Rickett, 1990, for a review), but we will concern ourselves in this paper with angular 
broadening alone. Following other authors in the study of the Galactic distribution of 
scattering material, we will make some simplifying assumptions. Apart from the fact 
that these assumptions are readily defendable, they also enable us to easily compare 
our results with others. 


2.1.1. The power law description 

Following Cordes, Weisberg and Boriakoff (1985), we describe fluctuations in the 
electron density by an isotropic power law: 

ftn.(q,*) = C* (*)r“ (1) 

which holds between lower and upper wavenumber cutoffs: qo < q < qi- C 2 (x) 
is the power in the fluctuation spectrum at some position along the line of sight 
(the changes in C 2 (x) take place on a much larger scale than the electron density 
fluctuations). There have been several studies specifically to determine a, qo and qt 
(Armstrong et al. 1981; Cordes et al. 1985, 1990; Gwinn et al. 1988; Spangler Sc 
Gwinn 1990). Determining these values is beyond the scope of our observations and 
so we use a = 11/3 throughout this paper. 

The power spectrum above will produce a phase structure term in a propagating 
wave through the medium. This will affect the measured visibility of a point source 
at a distance L, which is described by, 

V(p,L) = exp[-Dt(p)/2) (2) 

where p is the baseline length and D+(p) t the phase structure function, is given by 

D *(P) = 8r , r*A 2 jf j q[l - J 0 (qp)]P f „.{q,x)dqdx 

for the case of a plane wave (source at infinity). For a source embedded in the medium, 
it has been shown (Ishimaru, 1978, p. 418) that Jo(qp) is replaced by Jo(qpz/L) where 
(x/Z) is a dilution factor. For 2x/qi < p < 2v/qo, this results in, 

K(p,Z) = expHpM)° _J ] (3) 

where 


p e = ^A’r^ZMa)]- 1 '* 0 - 2 ) 


(4) 
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in which /(a) is an integration constant equal to 1.11835 for a = 11/3. S(L) takes 
into account the integration of C 3 (x) from source to observer. It equals 

S(L) = £ C' n (z)(l) a ~ 2 dz (5) 

when the observed source lies within the medium, or 

S(L) = ( L C 2 (z)dx (6) 

Jo 

when observing an extragalactic source. In this last case S(L) is by definition the 
scattering measure, SM. This quantity is directly usable in the case of an extra r 
galactic source. For a source embedded in the medium, as in this work, an additional 
assumption is needed before the SM can be used. If we make C 3 (z) constant, equa- 
tion (5) can be expressed as S(L) = (a — l)" l 5Af (Cordes et al., 1991a). If C 3 is not 
a constant then equation (5) should be used. Note that there is a strong weighting of 
scattering along the line of sight. The effect is that electron turbulence close to the 
source is weighted much less strongly [oc (z/L ) a ~ 2 ] than material more distant from 
the source. 


2.1.2. The very strong scattering regime 

The analysis embodied in equations (3) to (5) assumes that the relevant baselines (i.e. 
p ss p c ) are bracketed by the inner and outer scales of the electron density power 
spectrum: 

lx = 2irqi l <p<l 0 = 2 (7) 

In very strong scattering, the visibility function becomes narrower than the inner 
scale, i\. In this case, different expressions for the scattering measure hold. Other 
researchers (e.g. Spangler and Gwinn 1990; Moran et al. 1990; Molnar et al. 1990) 
have found evidence in favor of an inner scale l\ ~ 100 km, suggesting that any 
measured angular diameters 0 o b t smaller than about 100 mas at 1 GHz are influenced 
by Unite inner scale effects. 

There are several consequences of a finite inner scale which should be mentioned. 
First, very strong scattering that yields a very small spatial scale for the visibility 
does not require correspondingly small scale sizes in the scattering medium. Second, 
in very strong scattering, the visibility function is a Gaussian and the scattering 
diameter scales as 0 o b, oc A 3 . Thus 

V{p, L) = exp[— (/ ?/p e ) 2 ] (8) 

where it can be shown (e.g. Cordes and Lazio 1991) that 

Pe = [2*W e S(L) q }-°( 4 - a) -1 ] ~ 1/3 . (9) 

Third, the Iine-of-sight weighting is changed slightly, such that, for sources embedded 
in the medium, equation (5) is replaced with 

S(L) = £c 2 n (z)(j ; j t dz . 

The result of equation (6) still holds for extragalactic sources. 


( 10 ) 
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2.1.3. Resulting expressions 

In general we give the following result, expressed in observable quantities 


&obt 


S(L) = 


133 mas 




&obt 




5/3 


_ 128 mas ) 


11/3 

^GHi 


Oob. < Oc 


( 11 ) 


where 9„ 0 „ ss 0.16"(i/ G H.(fi/100 km)) -1 . 

To give the reader some idea about the numbers involved, we mention here that the 
values of C 3 (x) determined for the solar neighborhood are ~ 10 - 35 m~ 20 / 3 , resulting 
in an angular size of about 1 mas for an extragalactic source at 18 cm for a line of 
sight of 1 kpc through the medium. C 3 is typically 1 m ~ 30 ^ 3 in the inner Galaxy, 
resulting in up to 150 mas for objects seen through the inner Galaxy. 


2.1.4- Electron Densities and Emission Measures 


The electron density variations responsible for scattering also contribute to free-free 
radiation. The calculation of the emission measure EM involves n t = n e + 6n e with 
n 9 the local mean density and Sn 9 the local variation about the mean. To estimate the 
rms electron density Sn 9 one integrates Ps n ., the fluctuation spectrum. For a = 11/3 
and an outer scale to » l\, 


«"•(""*) = (&) " “ »° 3 ( w& s*) 1,1 (xU) 1/3 • (12) 


This will give a contribution to the EM, which may be written as 

2/3 


EMsh « SM = 10 274 pc cm -6 


SM 


a — 3 t0 F \lpcj V 1 kpc m -20 / 3 

(Cordes et al. 1991a, b). This corresponds to a minimum optical depth 

SM \ ( t 0 V' 3 g(u, T) 


) (13) 


r//SM - 10 3 54 (l kpc m - 2 °/ 3 )(lpc) 


1,2 

"gh* 


(14) 


where T 4 = (T/10 4 K) and g(u, T) = 1 + 0.085 In T 4 + 0.387 In ^qh*- Density fluctua- 
tions are bounded, locally, by 6n 9 < n 9 , so the actual EM and optical depth are at 
least a factor of two larger than the values in equation (13) and equation (14). 


2.2. oh/ir stars 

In order to fully appreciate why we can use OH/IR stars to study interstellar scatter- 
ing, we must first discuss the properties of OH masers in circumstellar shells. From the 
moment of their discovery it was suggested that the strong 1612 MHz spectral lines 
with their typical two-peaked signature are caused by maser emission in expanding 
circumstellar shells (Wilson he Barett, 1968). In the last two decades several aspects 
of these OH/IR stars have become clear. First, we now know that the underlying stars 
are objects on the Asymptotic Giant Branch. Having reached their final stages of evo- 
lution, these stars are losing a considerable percentage of their mass (< 10 -4 Mg/yr), 



resulting in a circumstellar shell which is gradually flowing away from the star. Dust 
condenses in this shell and is responsible for the prominent IR signature of these stars. 
Because of the large, cool, optically thick dust shell we see all the luminosity of the 
giant star being emitted longward of 5 pm. 

This IR emission causes a population inversion in the OH molecules, which are 
abundant in the shell. Provided that long coherent pathlengths are available, line 
emission in such a medium can develop into a maser. Precisely because the shell 
is expanding, the maser works most efficiently along the radial lines in the shell. 
Here all OH molecules in the shell can “see” one another; in more tangential lines 
the different parts of the shell cannot amplify each other’s emission because of a 
relative Doppler shift. Therefore every observer in any direction from the star sees 
the brightest emission from the extreme back and front of the shell, the outermost 
blue- and redshifted parts of the spectrum. 

This maser beaming, arising from the fact that the shell is expanding, results in a 
higher brightness temperature for the peaks of the spectrum than for the inner parts. 

From interferometric observations it is clear that the brightness temperatures of the 
inner parts of the spectrum are typically of the order of 10 4 K, while for the peaks 
of the spectrum this can be at least 10 8 K (see Alcock k Ross, 1988, and references 
therein). In several cases (Norris et al., 1984, Benson k Mutel, 1979, Bowers et al., 
1980, Reid et al., 1977) hotspots with angular sizes of the order of 10 mas are detected 
which have corresponding.10 9 - 10 10 K brightness temperatures. 


2.3. oh/ir stars as probes of interstellar scattering 

In this paper we interpret measurements of the angular sizes of the peak maser emis- 
sion from OH/IR stars as scatter broadening caused by interstellar scattering. Of 
course we can do this only when we can safely rule out that the intrinsic angular 
brightness distribution has no effect on our measurements, treating the underlying 
OH masers as point sources. For most sources in the survey we present here, we 
detect brightness temperatures of just 10 8 K for the peak of the spectrum, which we 
safely attribute to interstellar scattering. In this respect it is interesting to discuss 
the result by Diamond et al. (1988), who show that the angular sizes of the peaks in 
OH/IR stars increase with distance. This result indicates that interstellar scattering 
determines the sizes (and brightness temperatures) for OH masers at large distances, 
as had been suspected earlier (Burke et al., 1968, Boyd k Werner, 1972, Moran et al., 
1973, Bowers et al., 1980, Kent k Mutel, 1982, Gwinn et al., 1988, Kemball et al., 
1988). 

Because all these authors fail to And sources with angular sizes less than 20 mas 
more distant than 1 kpc, there is good reason to believe that the maximum bright- 
ness temperature measurements for OH masers are frequently limited by interstellar 
scattering. With this in mind, we feel safe to use all measurements of the angular size 
above 50 mas as a measurement of interstellar scattering for sources at the distance 
to the Galactic center. This might also be valid for sources with smaller sizes, but to 
be conservative these should be interpreted as upper limits on angular broadening. 


TABLE 1 

List Op Observed OH/IR Stars 


Source 

Name 

(1) 

RA(1950) 

(2) 

DEC (1950) 
(3) 

Peak Velocities 
(km s -1 ) (km a -1 
(4) 

Peak Fluxes 

)(Jy) (Jy) 

(5) 

IRAS 

Designation 

(«) 

OH 357.849+9.744 ... 

... 17 00 49.34 

-25 12 23.3 

-31.6 

-11.2 

1.7 

0.5 

17008-2512 

OH 0.125+5.111 ... 

.. 17 23 22.90 

-26 02 24.1 

-151.4 

-135.6 

2.3 

1.1 

17233-2602 

OH 353.298- 1.537 ... 

. . 17 31 36.51 

-35 23 55.9 

48.5 

74.6 

3.1 

3.6 

17316-3523 

OH 359.140+1.137 ... 

. . 17 35 56.95 

-29 02 24.7 

-145.7 

-127.5 

2.4 

1.8 

17359-2902 

OH 359.564+ 1.287 ... 

. . 17 36 24.70 

-28 36 10.0 

17.9 

45.2 

1.7 

0.5 

17364-2836* 

OH 355.641-1.742 ... 

. . 17 38 34.34 

-33 32 13.9 

-245.7 

-224.2 

3.3 

5.1 

17385-3332 

OH 355.897- 1.754 ... 

..17 39 16.11 

-33 19 36.1 

-58.0 

-30.8 

1.1 

1.4 

17392-3319 

OH 0.892+1.342 ... 

.. 17 39 24.74 

-27 27 02.0 

-123.3 

-93.7 

49 

0.8 

17393-2727 

OH 1.212+1.258 ... 

.. 17 40 29.90 

-27 13 25.3 

30.4 

57.7 

5.8 

4.9 

17404-2713 

OH 359.517+0.001 ... 

. . 17 41 16.555 

-29 19 35.19 

-145.6 

-111.5 

0.9 

1.5 


OH 359.762+0. 120 ... 

. . 17 41 24.132 

-29 03 20.94 

-21.0 

9.6 

1.8 

3.3 


OH 1.369+1.003 . . . 

. . 17 41 50.64 

-27 13 30.4 

-27.4 

3.3 

8.7 

6.8 

17418-2713 

OH 359.581 -0.240 ... 

. . 17 42 22.155 

-29 23 56.71 

-103.6 

-75.2 

3.1 

0.9 


OH 359.880-0.087 ... 

. . 17 42 29.643 

-29 03 53.70 

-34.6 

-14.2 

6.6 

2.0 


OH 359.986-0.061 ... 

..17 42 38.757 

-28 57 40.35 

-7.4 

50.5 

1.2 

1.0 


OH 0.319-0.040 ... 

. . 17 43 21.902 

-28 39 58.32 

57.3 

92.2 

1.3 

3.4 


OH 0.814+0.179 ... 

. . 17 43 41.639 

-28 07 47.81 

56.2 

88.8 

0.5 

1.9 


OH 0.334— 181 ... 

.. 17 43 56.62 

-28 43 39.1 

-356.4 

-326.9 

3.5 

2.9 


OH 5.026+1.491 ... 

. . 17 48 26.09 

-23 50 56.3 

106.7 

135.1 

3.0 

1.2 

17484-2350* 

OH 3.234-2.404 ... 

.'. 17 59 18.87 

-27 21 58.8 

-45.2 

-20.2 

4.6 

0.2 

17593— 2721* 


a . 17363-2833 b . 17486-2345 C. 17596-2716 


3. Observations 
3.1 SOURCE SELECTION 

In order to map the distribution of interstellar scattering about SgrA*, we require 
a large sample of bright OH/IR stars. Fortunately, there is a strong concentration 
of OH/IR stars within one degree of the Galactic center (Lindqvist et al., 1991). 
Approximately half of the 20 OH/IR stars used in our work (see table 1) were taken 
from the Lindqvist et al. (1991) survey in which 134 sources were discovered in the 
central 1.5 square degrees of our Galaxy. Since the VLA 1 was used for their survey, 
subarcsecond positions were already available for our work. For sources that have 
|6| > 0.5°, the all-sky survey of te Lintel Hekkert et al. (1991) was used instead. 
The best positional information in the te Lintel Hekkert et al. (1991) survey was 
based upon a match between the OH/IR star (detected at 1612 MHz) and its nearest 
IRAS point source with the proper colors. These positions were only accurate to ±15 
arc sec and so a short VLA astrometry program had to be carried out before the VLBI 
observations could be done. 


1 The Very Large Array (VLA) and the Very Long Baseline Array (VLBA) are operated by the 
National Radio Astronomy Observatory under cooperative agreement with the National Science 
Foundation. 
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3.2. ASTROMETRY 

On June 21, 1990 the VLA was used in the hybrid A/B array configuration and two 
polarizations of 64 spectral line channels each were recorded across a bandwidth of 
391 kHz (or 6.1 kHz per channel). The data were reduced within the AIPS package 
following standard practice. Table 1 contains the results of this work along with 
the positions of the lower latitude sources from the Lindqvist et al. (1991) catalog. 
Column (1) lists the source name given with a preface of “OH” , followed by its galactic 
coordinates, while columns (2) and (3) are the right ascension and declination of each 
source, respectively. Column (4) gives the velocities for the blue- and red-shifted peaks 
of the double-peaked profiles. Similarly, the peak flux densities for these components 
are given in column (5). Because OH/IR stars are known to be highly variable these 
values should only be taken as approximate. 

We observed 15 OH/IR stars from the catalog by te Lintel Hekkert et al. (1991); 
13 were detected. 17338—2744 and 17392—2722 were not detected. For all but one of 
the 13 detected sources there is an IRAS source name in column (6). OH 0.334—0.181 
was paired with 17439—2845 by te Lintel Hekkert et al. (1991), but at our new VLA 
position, 1.5 arcminutes away, there is no known IRAS source. This is likely due to 
confusion in the IRAS data, because of the low latitude. Another 3 OH/IR stars out 
of these 13 were found to be associated with an IRAS point source other than the 
one that te Lintel H>.kkert et al. (1991) assigned. The original IRAS names of these 
sources are found just below Table 1. Apart from these 6 cases we believe the other 
IRAS assignments to be correct. The poor detection statistics should not be taken 
as typical for the te Lintel Hekkert et al. (1991) survey, because the stellar density 
is very high in this region. There are no IRAS designations for the Lindqvist et al. 
(1991) sources in table 1, owing to the severe confusion that occurs in the infrared at 
these lower latitudes. 


3.3. THE DISTANCE PROBLEM 

The selection of appropriate OH/IR stars from these catalogs was governed primarily 
by considerations of sensitivity. In order to make a detection within the coherence 
time on a VLA-VLBA baseline with the characteristically small linewidths of these 
masers (~ 6 kHz), we must select stars with a peak flux above 1.5 Jy. This condition 
was relaxed somewhat in a few instances along interesting lines of sight. In addition, 
as a check, we included one of the sources from van Langevelde k Diamond (1991). 

To simplify any comparison of the scattering properties between these different lines 
of sight, it is important to choose the sample of OH/IR stars so that they are located 
at or beyond the Galactic center. From their flattened distribution and kinematics we 
know, at least in a statistical sense, that the Lindqvist et al. (1991) survey contains 
stars at the Galactic center. For the te Lintel Hekkert et al. (1991) sample, the 
situation is a bit less certain because the density contrast towards the center of the 
Galaxy is less pronounced than in the Lindqvist et al. (1991) sample. Fortunately, the 
IRAS data yield a useful discriminant. While the 1612 MHz maser emission can vary 
widely from one star to the next (Gaylard et al., 1989, Rottgering, 1989, Dickinson, 
1987), the variation in bolometric stellar luminosity is relatively more well-behaved. 
This enables us to obtain some distance information about individual objects. 

Thus, obvious foreground stars can be eliminated from our sample by choosing only 
those stars whose IRAS fluxes are weak. We have estimated a distance from the IRAS 
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12 and 25 /im bands using a standard luminosity of 15000 L©, which is a reasonable 
estimate of the luminosity of these OH emitting AGB stars (e.g. van Langevelde et 
al, 1990). For this estimate we have used the bolometric correction for OH/IR stars 
given by van der Veen ic Breukers (1989). We eliminated all sources whose distances 
calculated in this manner were less than 8 kpc. Work to derive the IR fluxes for 
the low latitude sources in Lindqvist et al. (1991) is in progress (Blommaert, private 
communication) . 


3.4. VLBI OBSERVATIONS 

The observations were made on two separate occasions, October 18, 1990 and February 
1, 1991. The phased VLA and all the working antennas of the VLBA were used in both 
cases. This included antennas at Pie Town (PT), Los Alamos (LA), Kitt Peak (KP), 
and Fort Davis (FD). The point source calibrator 1748—253 was observed for several 
minutes each half hour to re-phase the elements of the VLA, and the bright source 
NRAO 530 was observed every 4 hours to act as a “fringe finder” during correlation. 
Each OH/IR star was observed for 10-12 minutes at a time, with two scans made at 
different hour angles so as to sample a full range of baseline lengths. 

The data were recorded on Mkll tapes with a total bandwidth of 500 kHz. The 
tapes were processed on the NRAO Mkll correlator, splitting the bandwidth into 96 
spectral line channels of 5.2 kHz each (0.97 km s -1 ). After correlation the data were 
further processed in AIPS following procedures outlined by Diamond (1989). The 
bandpass correction was made using the total power spectrum from each antenna for 
NRAO 530 and the data were corrected for Doppler tracking. In the first session, the 
amplitude calibration was done in the same manner as explained in van Langevelde & 
Diamond (1991) by using the monitored system temperatures and estimated antenna 
gain factors, in much the same way as is done with continuum VLBI. In the second 
session, however, the VLA data showed signs of severe broadband contamination, 
most likely from the GLONASS network of satellites. A full amplitude calibration 
was impossible. Instead we either estimated calibration constants from a priori infor- 
mation, or used only VLA-VLBA baselines for further analysis. Because all VLBA 
antennas are essentially identical, we expected VLA-VLBA baselines to be equally 
sensitive within 10%, as was found during the first run. 

A fringe fit was done for our calibrator, NRAO 530. We also had observations of 
1748—253 every half hour to re-phase the VLA, so we were able to use it to fringe fit. In 
this way we obtained reasonable solutions for delay and fringe rates for all antennas 
during the entire observation period. We then continued with a fit of the residual 
fringe rate on the OH maser sources, setting a window on the spectral channels in 
which we expected the peak fluxes. We set the threshold for an acceptable detection 
to 4<r, where we had limited the search for a fringe to the spectral channel where we 
expected a signal. By performing a residual fringe fit we gave up the absolute position 
information on the OH source. 

In principle, the final step in obtaining the results we are interested in is straight- 
forward. We have the spectrum of the source from a number of baselines and a “zero 
spacing” measurement from the autocorrelation spectrum on any of the telescopes. 
The autocorrelation spectrum on the VLA is most often used for the latter, because it 
has the highest signal-to-noise and benefits from the lower confusion provided by its 
small synthesized beam. At this point we simply fit a Gaussian profile to the visibility 
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data of the brightest channel. We determine 6 0 i,, the full width at half maximum of 
the corresponding Gaussian in the map plane, directly from this value. 

Our basic dataset and fitting procedure is illustrated in figure 1 where we show the 
results for OH 0.125+5.111. Although this is a relatively weak source, it was detected 
on all the baselines present (PT-LA is missing, because of a technical problem). Note 
that this particular source is not heavily scattered, while for most of the sources in 
our sample, detections were made on only the shortest and most sensitive baselines. 



Figure 1. The amplitude of OH 0.125+5.111 a* • function of baaeline. The top row thowz the spectra 
a a detected on baselines of increasing length towards the right. The bottom shows the amplitude of 
the brightest, blue-shifted peak as a function of baseline length together with the fitted model far 
a spherical gaussian. Because we ^t"inr a single component and the data have been subjected to 
fringe fitting we d is re ga r d the phase information and take it to be zero on average. 


In table 2 we list the measured angular sizes of all sources which were detected 
on one or more baselines. Also given in table 2 is the log of the scattering measure 
for a scattering medium assumed to have constant C^(z). The last column in table 2 
is the angular distance from SgrA*. Seven of our 20 sources were not detected on any 
of the measured baselines but were seen in the autocorrelation spectra. For three of 
these sources we did not have data for our shortest baselines, because of a failure at 
PT, so it is perhaps not surprising that we did not detect these. 

We do not give errors on the measurements for each source in table 2. An accurate 
estimate of the uncertainties in the angular sizes is difficult to give, especially for the 
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larger scattering sizes. They can generally not be estimated from the formal errors in 
the fit because the luge angular sizes are often based on detections of a single baseline 
and are therefore a strong function of calibration and the detection threshold. Taking 
this into consideration we estimate the errors on the large scatter sizes to be as high 
as 30%. For the small scatter sizes we think we can rely upon the formal errors in the 
fit being on the order of 5% - 10%. 

TABLE 2 

Scattering Diameter Measurements 
Of OH/nt Stars At 1612 MHz 

Source 0 O », logSM* ASgrA* 

Name (mas) (kpc (arcmin) 


(1) (2) (3) (4) 


OH 353.298- 1.537.. 

... 360 

2.17 

408.6 

OH 355.641 -1.742.. 

... 329 

2.09. 

277.5 

OH 355.897—1.754 .. 

... 348 

2.14 

263.5 

OH 357.849+9.744.. 

... 157 

1.45 

600.6 

OH359.140+1.137 . . 

... 269 

1.92 

85.8 

OH 359.564+1.287.. 

... 81 

0.86 

83.2 

OH 359. 762+0. 120.. 

... 495 

2.45 

14.8 

OH 0.125+5.111 .. 

... 35 

0.25 

309.6 

OH 0.892+1.342 .. 

... 57 

0.60 

100.9 

OH 1.212+1.258 .. 

... 258 

1.88 

109.1 

OH 1.369+1.003 .. 

... 112 

1.16 

106.2 

OH 3.234-2.404 .. 

... 141 

1.36 

242.8 

OH 5.026+1.491 .. 

... 68 

0.73 

318.5 


f Assuming a medium with uniform Cj and f j = 100 km. 


3.5. KNOWN LINES OF SIGHT 

To aid in our discussion we can enlarge the number of sampled lines of sight by 
including measured angular diameters from the literature. The 4 OH/IR stars of 
van Langevelde tc Diamond (1991) were taken from the Lindqvist et al. (1991) sample 
and were observed with exactly the same telescope array and procedures as described 
above. We also use the measurements of OH and HjO masers in the SgrB2 complex 
(Gwinn et al. 1988, Moran, 1968). In addition, two radio transients have been 
discovered close to the Galactic center, both of which have had angular broadening 
measurements made over a range of wavelengths (Zhao et al. 1991, Davies et al. 
1976). Finally, of course, we include SgrA* itself (Lo et al., 1985). For these last 
four sources it is known that their angular sizes scale quadratically with wavelength, 
consistent with their sizes being determined by interstellar scattering. The authors 
in each of the above-mentioned references also argue that their sources are located at 
the Galactic center. These extra eight lines of sight are presented in table 3 in much 
the same format as table 2. 
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TABLE 3 

Additional Angular Broadening Measurements Scaled To 1612 MBs 


Source 1 b 9 logSM ASgrA* Ref. 

Name (degrees) (degress) (mas) (kpc m~ 10 / 3 ) (arcmin) 

(1) ffl (3) (4) (5) (8) (7) 

OH 359.762+0.120 359.762 0.120 531 2.51 14.8 1 

OH 359.938-0.077 359.938 -0.077 571 2.57 1.9 1 

GCT 359.939 -0.055 500 2.46 0.6 2 

SgrA* 359.944 -0.046 520 2.49 0.0 3 

A 1742-28 359.930 -0.043 530 2.51 0.9 4 

OH 359.954-0.041 359.954 -0.041 591 2.60 0.7 1 

OH 0.190+0.036 0.190 0.036 103 1.08 15.6 1 

SgrB2 0.7 0.01 96 0.98 45.5 5 


(1) van Langevelde and Diamond (1991) (2) Zhao cl al (1991) (3) Lo cl at. (1985) (4) Davies 
cl al (1978) (5) Moran (1968). 

4. Discussion 

OH 359.762+0.120 was observed both here and by van Langevelde k Diamond (1991), 
and the 0 o i, agree within the errors. 

Figure 2 shows the distribution of scattering sizes in galactic coordinates for the 
data taken from Tables 2 and 3. The sizes of the open circles are directly proportional 
to 9 oh ,. The most striking aspect of this figure is the concentration of large scattering 
sizes near SgrA*. This is more evident in Figure 3 where we show the observed 
scatter size as a function of distance from SgrA*. There are 7 measurements within 
±15 arcmin of SgrA* with an average 0 o t, of 525 mas. In addition, we draw attention 
to the lower limits on 0 O ». (upper limits on correlated flux) found on the left hand 
side of figure 3. These are the 7 sources for which we failed to detect VLBI fringes. 
They are not randomly distributed about the range of sampled distances from SgrA*, 
as would be expected if our non-detections were due to errors in pointing or poor 
sensitivity. It suggests that these lines of sight are more scattered still and that they 
are all but resolved by even our shortest baselines. 

There are prominent variations in in figure 2 over short angular scales. Beyond 
the central region of 15 arcmin radius, 9 0 tt drops precipitously to a value of 100 
mas. Two sources are responsible for the drop, SgrB2 and OH 0.190+0.036. In the 
past this result for SgrB2 has been used to define the outer limit for the size of 
the scattering region around SgrA* to be 45 arcmin, or a maximum length scale 
of 150 pc (Backer, 1978, Ozernoi k Shishov, 1977). With OH 0.190+0.036 we can 
constrain the outer extent further still to 15 arcmin, or a radius of 40 pc. However, 
we caution that there are no IRAS colors for this star; thus the possibility exists 
that it is a foreground object. At 6 ~ 1.5 degrees there are two additional examples 
of prominent variations taking place on short angular scales. Two pairs of sources, 
OH 359.140+1.137 and OH 359.564+1.287, and OH 1.212+1.258 and OH 0.892+1.342, 
are believed to be located at or beyond the galactic center on the basis of their IRAS 
fluxes. In each pair the ratio of 9„ht varies by a factor of 3-5 over angular scales ot 20-25 
arcmin, which corresponds to order-of-magnitude variations in SM, over a maximum 
transverse distance of only 50-65 pc (but an unknown distance along the line of sight). 
We also note that the data do not prove that scattering decreases monotonically from 
the Galactic center, because some of the lower bounds are more than a factor of two 
larger than actual detections for similar angular distances from the Galactic center. 
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Figure 2. Scattering sizes for all the lines of sight near the Galactic center. The sizes of the circles 
indicate the amount of observed angular broadening. 


4.1. COMPARISON WITH OTHER SCATTERING MEASUREMENTS 

To put the Galactic center source measurements into context, we show in Figure 4 
their scattering measures along with those of pulsars, other maser sources, and one 
active galactic nucleus (this figure is copied from Cordes et al. 1991a) for the area 
—30° < l < 30®, —10® < b < 10®. As can be seen, SM increases with distance much 
faster than would be expected if Cj| were constant and equal to the value inferred for 
the local interstellar medium. The rise occurs within the first few kpc from the Sun. 
The solid line in Figure 4 displays SM predicted from a best fit model (Cordes et 
al. 1991a) to more than 200 lines of sight. The model consists of an outer Galaxy 
component that is nearly constant in galactocentric radius and has a z scale height of 
1 kpc; an inner Galaxy annulus is centered on R as 4 kpc, has a e -1 width of 2 kpc, 
and a scale height of 100 pc. In terms of C*, the inner Galaxy component is about 
2000 times stronger than the outer component. 

The model approximates SM data in other directions (e.g. toward the galactic poles 
and toward the galactic anti-center) with relatively little scatter. The large amount 
of scatter in Figure 4 is due to the small scale height (100 pc) of the inner Galaxy 
component, corresponding to a galactic latitude of only 0.7® at 8.5 kpc distance. The 
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Figure 3. Plot of scatter sizes with distance from SgrA* 

dotted line shows the model’s SM for jf| = 15°, 6 = 0° while the dashed line shows the 
model looking out of the galactic plane at |6| = 5°. Thus, points below the solid line 
may be understood as showing smaller SM because of the latitude dependence built 
into the model, while points above the line correspond to lines of sight with enhanced 
scattering (i.e. scattering in excess of that predicted by the model). For example, 
the masers inside NGC6334 show enhanced scattering (as indicated) and the AGN 
behind it also shows enhanced scattering. The difference in SM for the background 
AGN and the maser is caused by the same leverage effect presented in equations (4) 
and (5). 

The most important feature of Figure 4 for the discussion here is the obvious en- 
hancement of points for the sources with attributed distances of 8.5 kpc. The figure 
shows clearly that the scattering measures of the galactic center sources are anomalous 
with respect to other lines of sight toward the inner Galaxy. We note, however, that 
the appearance of Figure 4, namely the sudden rise oi SM from that of the model 
to about 10 3 4 kpc m~ 30 ^ 3 , does not prove that the enhanced scattering occurs at 
the distance of SgrA*; none of the other sources plotted in the figure is at the same 
longitude and latitude, so the onset of the large enhancement could in fact occur, at 
a distance smaller than 8.5 kpc. 

The line of sight towards SgrA* and the 30 arcmin diameter region around it does 
show an enhanced level of scattering, even when compared to lines of sight in the 
immediate vicinity. Where, then, does this turbulent medium originate and what can 
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Distance [kpc] 


Figure 4. Plot of scattering measure against distance from the sun. The plotted points represent 
data presented in Tables 2 and 3 and data compiled by Cordes et al. (1991a) that include pulsar 
measurements and data on scattered AGNs from Fey et al. (1991). AGNs are plotted arbitrarily 
at a ‘distance’ of 30 kpc. The solid line represents SAf calculated from the model of Cordes et aL 
for / = 0, i = 0; the dotted line corresponds to jf| = 15°, 6 = 0; the dashed line corresponds to 
< = 0,161 = 5°. 


we say about its location relative to SgrA*? We consider two possible configurations 
to account for the turbulence: an anomalous scattering region close to SgrA* and an 
unrelated turbulent region distributed in the interstellar medium towards SgrA*, but 
not necessarily at the same distance as SgrA*. 

In the discussion we will adopt the formulas derived for very strong scattering, 
although the observations presented here do not yield any direct information on the 
presence of an inner scale. However, as pointed out in the introduction, it seems 
likely that observations with 9 0 bt > 100 mas at 1 GHz are affected by inner scale 
effects. Furthermore, the observed angular size of SgrA* <x A 3 seems to support this 
conclusion. 


4.2. ANOMALOUS SCATTERING LOCAL TO THE GALACTIC CENTER 

Ozernoi & Shishov (1977) were the first to suggest that the scattering seen for SgrA* 
occurred local to the source. They calculated C 3 for a single scattering screen located 
at several different positions along the line of sight, and then, assuming an outer scale 
to, they computed Sn e . They concluded that the ratio Sn t /n e was “reasonable” only 
when the scattering region was located 10 pc from SgrA*. Backer (1978) indepen- 
dently concluded that the scattering towards SgrA* was anomalous and gave the same 
distance of 10 pc for the location of the scattering screen. 
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Our observations also may be interpreted in terms of a local origin for the scattering. 
As mentioned above, the zone of enhanced scattering appears to have a well-defined 
angular extent of 30 arcmin, approximately centered on SgrA*, and there are no other 
lines of sight in our sample with comparable 9 0 t» ■ There are strong phenomenological 
reasons for believing in a local origin, as the Galactic center is arguably one of the 
more energetic environments in our Galaxy. 

We develop a simple model to examine this hypothesis, and use it to determine 
the physical properties of the scattering regions. In our model we assume there is 
an extended medium with C 3 (r) = C 3 0 , and a region of enhanced scattering with 
C 2 (*) = C 2 oc which is distributed in a sphere centered on SgrA* with a radius Rqc- 
This model is fitted to the data and the best-fit parameters are determined using the 
X 3 method. A distance of 8.5 kpc is assumed for the distance of each source. This 
results in a fit of C 3 >0 = 6.3 m -30 / 3 , C% ac = 2.1 • 10 8 m -2 °/ 3 and Rqc = 54 pc. This 
value for C 2 GC differs from that estimated by van Langevelde It Diamond (1991) for 
a similar geometry, because here we include inner scale effects. 

It is tempting to try to link the enhanced scattering to a known region at the 
Galactic center. In this respect we note that the size of 50 pc coincides with the size 
of the “magnetosphere” of the Galactic center (Heyvaerts, 1989). 


4.3. ANOMALOUS SCATTERING DUE TO CLUMPS AWAY FROM THE GALACTIC CENTER 

It is well known that some lines of sight show enormously more scattering than neigh- 
boring ones, such as the Vela pulsar and Cygnus X-3. Moreover, the HII complex 
NGC6334, at t — 351°, 6 = 0.7°, and D a 1.7 kpc, causes scattering of a background 
extragalactic source that is larger than the scattering of the Galactic center sources 
discussed here (Moran et al. 1990). Cordes et al. (1985) discussed a model for the 
scattering medium which has a smoothly distributed component along with clumps of 
very strong turbulence. Clumps of parsec size were inferred to have a volume filling 
factor of about 10~ 4 , corresponding to a mean free path for intersection of a clump of 
about 5 to 10 kpc. Using a more extensive data set, a model with two smooth compo- 
nents has been devised (Cordes et al. 1991; as shown in Figure 4) that also must be 
augmented with clumps of turbulence with the same properties as the earlier model. 
Because the mean free path for intersecting a strong clump is less than or about equal 
to the distance to the Galactic center, it remains plausible that a similar clump may 
account for the enhanced scattering of the OH/IR masers and other sources in the 
Galactic center. 

Previous authors have reported variations of scattering on angular scales of a few 
degrees or more (Fey et al. 1989, Alurkar et al. 1986, Dennison, et al. 1984). With our 
much finer grid of sampled points we are able to put an important limit on the typical 
size of a scattering clump, namely w 50 pc • (D/8.5 kpc) where D is the distance of 
the clump from the Sun. The upper bound of 50 pc is similar to clump sizes discussed 
by Cordes, Weisberg, and Boriakoff (1985) and identical to that used by Fey et al. 
(1991) in their detailed model simulations. 


4.4. IMPLICATIONS AND PREDICTIONS 
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4-4-1- Extragalactic sources 

The model of enhanced scattering local to the Galactic center gives an interesting 
prediction for the observations of compact extragalactic sources, particularly those 
with flat spectra. Putting such a source behind the scatterer, within 15 arcminutes 
of SgrA*, we expect a size of 213" at 1.6 GHz. For a scattering medium closer to 
the earth, this prediction comes out smaller, depending on the exact location of the 
scatterer. The best approach would be to make observations at higher frequencies (e.g. 
near 20 GHz with the VLA) where the angular broadening is easier to measure (if the 
scattering arises in the vicinity of the Galactic center). Note that angular broadening 
decreases faster with an increase in frequency than does the angular resolution of a 
telescope. 


4-4-2- Physical Implications 

The location of the region that causes enhanced scattering of the Galactic center 
sources has a large influence on the amount of free-free absorption to be expected. 
Consider a clump model for the variation of C% with distance from the Galactic center: 


Cl(x) = 


Clo + Cl 


n,GC 


U B, 0 


\z-L ac \<AL/2 


otherwise 


(15) 


where L ge is the distance of the clump from the galactic center and A L is its thickness. 
Here, for simplicity, the value of C* outside the clump, C* >0 , is assumed to be constant 
and much less than that inside the clump, C* QC . It may be shown, for A L/L <. 1, 
that 


S(L) a -5 Mq + SAf ac 


(¥)'• 


The total SM is 


SM = SMo + SMqc ** 


(£)' 


S(L). 


(16) 


(17) 


Introducing S(L) from equation (17) into equation (14) and using equation (11) we 
obtain 


T 11 


~ (jl v ( y ( l y /3 ( jo_\ 

~\LgcJ \7831.6ma s/ \100km/ \lpc/ 


-2/3 




9(v,T) 


'GB* 


r 4 15 


(18) 


where vs t GH$ is the frequency at which the angular broadening measurement is 
obtained and i/qh» is the frequency at which the optical depth is evaluated; l\ is the 
inner scale, while is the outer scale. Using as 500 mas at ug = 1.6 GHz we 
obtain 


te) 


~ 2/3 ^gh, 9(v,T) 


Tl s 


( 19 ) 


To have r// < 1 at 1.6 GHz requires that 


Lgc 

L 





It is possible to adjust the free parameters, particularly the temperature, to reduce 
the coefficient in equation (20) from 0.10 to 0.01. It is also arguable that optical depth 
unity occurs well below 1.6 GHz, because the spectrum for Sgr A* shows no break 
above 1 GHz (Lo et al. 1985). However, there is no question that the optical depth 
is greater than unity at 327 MHz, because Sgr A* is absent in VLA images at that 
frequency (Anantharamaiah et al. 1991). Indeed, Pedlar et al. (1989) conclude that 
Tjj in front of Sgr A* is > 2.5 (at 327 MHz), though it cannot be significantly larger 
than this value, otherwise the spectrum would turn over above 1 GHz. This would 
imply that 


< 0.32 

Lt 



( 21 ) 


Therefore, observations of Sgr A* suggest that, if the scattering and absorbing material 
are one and the same, then (1) the outer and inner scales of the scattering medium 
must be nearly equal to the fiducial values of 1 pc and 100 km, respectively; and (2) 
that the scattering medium is situated quite far (> 0.85 kpc) from the Galactic center. 


4-4.8. Implications for the variability of Sgr A* 

The observations presented here give some idea about the electron density fluctuations 
close to Sgr A*. Even in the model in which the scattering medium is located in the 
vicinity of the Galactic center, we find, from the distribution of scattered sizes, that 
it extends to 50 pc. A scatterer local to the Galactic center is also assumed in the 
discussion by Zhao et al. (1988). They infer that refractive scattering can explain the 
variability of SgrA* if the medium is moving with a velocity >1000 km s'" 1 relative to 
the source. Our observations make it unlikely that the diffractive scattering described 
in our paper and the refractive scattering postulated by Zhao et al. (1988) are caused 
by the same medium, because at a distance of 50 pc from Sgr A* the velocities are 
reasonably well known and significantly below 1000 km s'" 1 . 


4 . 4 . 4 . Possible associated magnetic irregularities 

Though the origin of electron density variations is unknown, it is plausible that they 
are related to, if not caused by, turbulence in the magnetic field such that 

£•<(¥)" 

with (al and n = 2 for MHD turbulence and n = 1 for obliquely propagating MHD 
waves (Cordes et al., 1990). If 6n e /n e is of order unity, as it seems to be for lines of 
sight through the Cygnus region and for other galactic plane lines of sight (e.g. Lazio 
et al. 1990, Clegg et al. 1991), then the energy density in magnetic field irregularities 
(smaller than about 1 pc) is equal to the energy density in the large scale field. 
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5. Conclusions 

By simple detection experiments of OH masers around the Galactic center we show 
that lines of sight to the Galactic center are subject to severe interstellar scattering. 
We also show from changes in the scattering taking place on angular scales of #s 20' 
that the upper limit on the size of a typical scattering clump is approximately 50 pc. 

We consider two alternate interpretations for the location of the enhanced turbu- 
lence. Either it arises local to the Galactic Center, and is directly related to the violent 
processes going on within or it is due to a clump of enhanced turbulence along the line 
of sight, and has no physical relationship to the Galactic Center. Our data can be fit 
to either model equally well but the former model seems to conflict with observations 
of free-free absorption towards SgrA, in which the same free electrons should play a 
role in scattering. Further observations are suggested to help resolve this controversy. 
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